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Cancellation of infrared divergences at NNLO 
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Abstract. Perturbative calculations at next-to-next-to-leading order for multi-particle final states require 
a method to cancel infrared singularities. I discuss how to setup the subtraction method at NNLO. 

PACS. 12.38.Bx Perturbative calculations 



1 Introduction 

The next generation of collider experiments will hunt for 

, the Higgs and other yet-to-be-discovered particles with 
increased luminosity and experimental precision. The in- 

[ creased experimental precision has to be matched by an 
improvement in the accuracy of theoretical predictions. 
Theoretical predictions are calculated as a power expan- 
sion in the coupling. Higher precision is reached by includ- 
ing the next higher term in the perturbative expansion. 
The experimental needs are numerical programs which 
yield predictions for a wide range of observablcs. Urgently 
needed are therefore fully differential next-to-next-to-lead- 

, ing order (NNLO) programs. Compared to certain spe- 
cific NNLO prediction for inclusive observables, these pro- 
grams are flexible and allow to take into account compli- 

■ cated detector geometries and jet definitions. The only re- 
quirement on the observable is infrared-safety. At NNLO 
this implies that whenever a n + 1 parton configuration 
Pi,...,p n +i becomes kinematically degenerate with a n par- 
ton configuration p' x ,...jp' n we must have 

O n+1 (p 1 ,...,p n+1 ) -> n (pi,...,p^). 

In addition, we must have in the double unresolved case 
(e.g. when an+2 parton configuration p\,...,p n +2 becomes 
kinematically degenerate with a n parton configuration 

Pl,-,Pn) 

O n+2 {pi,---,Pn+2) -> O n (Pi,...,j4). 

To construct such NNLO programs the following ingredi- 
ents are needed: 

- The scattering amplitudes. This implies in particular for 
a NNLO program the calculation of the relevant two- 
loop amplitudes. There has been substantial progress 
in this field in the past years. The state-of-the-art is 
that all two-loop-amplitudes, which are needed most 
urgently, are now known [1-13]. 
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- A NNLO program requires a method to cancel infrared 

divergences. Loop amplitudes, calculated in dimension- 
al regularization, have explicit poles in the dimensional 
regularization parameter e = 2 — D/2, arising from 
infrared singularities. These poles cancel with similar 
poles arising from amplitudes with additional partons, 
when integrated over phase space regions where two 
(or more) partons become "close" to each other. How- 
ever, the cancellation occurs only after the integration 
over the unresolved phase space has been performed 
and prevents thus a naive Monte Carlo approach for a 
fully exclusive calculation. It is therefore necessary to 
cancel first analytically all infrared divergences and to 
use Monte Carlo methods only after this step has been 
performed. 

- The final numerical computer program, which evalu- 
ates the remaining phase space integrals, requires sta- 
ble and efficient Monte Carlo methods for this integra- 
tion. 

In this talk I focus on the cancellation of infrared diver- 
gences [14,15]. In the next section I review general meth- 
ods at NLO. In sect. 3 I discuss the subtraction method 
at NNLO. Sect. 4 is devoted to one-loop amplitudes with 
one unresolved parton. 



2 A review of the subtraction method at NLO 

Infrared divergences occur already at next-to-leading or- 
der. As a simple example two diagrams contributing to 
the NLO corrections to e + e~ — > 2 jets are shown in fig. 1. 
The diagrams are divided into virtual and real corrections. 
The virtual corrections contain the loop integrals and can 
have, in addition to ultraviolet divergences, infrared diver- 
gences. For one-loop amplitudes the IR divergences man- 
ifest themselves as explicit poles in e up to 1/e 2 . For each 
IR divergence in the virtual corrections there is a cor- 
responding divergence with the opposite sign in the real 
emission amplitude, obtained from the integration over 
the phase space region where some particles become soft 
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Fig. 1. Cancellation of divergences between virtual and real 
corrections at NLO. 

or collinear (e.g. unresolved). In general, the Kinoshita- 
Lee-Nauenberg theorem guarantees that any infrared-safe 
observable, when summed over all states degenerate ac- 
cording to some resolution criteria, will be finite. However, 
the two contributions (virtual and real) live on different 
phase spaces and prevent a naive Monte Carlo approach. 
At NLO, general methods to circumvent this problem are 
known. This is possible due to the universality of the sin- 
gular behaviour of the amplitudes in soft and collinear 
limits. Examples are the phase-space slicing method [16- 
18] and the subtraction method [19-23]. I briefly review 
the subtraction method here. The NLO cross section is 
given as the sum of the virtual and real corrections: 

a NLO = J da R + J da v . 

n+l n 

If one can find an approximation term du A such that 

— da A has the same point-wise singular behaviour in D 
dimensions as da R itself, 

— da A can be integrated analytically in D dimensions 
over the one-parton subspace leading to soft and col- 
linear divergences, 

then one can add and subtract this term as follows: 




Since by definition do has the same singular behaviour 
as da R , da A acts as a local counter-term and the com- 
bination [da R — da A ) is integrable and can be evaluated 
numerically. Secondly, the analytic integration of da A over 
the one-parton subspace will yield the explicit poles in e 
needed to cancel the corresponding poles in da v . 

3 The subtraction method at NNLO 

The following terms contribute at NNLO: 

da {1) -(a {0} * A (1) +A (1) * A (Q) )dd>M, 
d&h ' — ^ A„ ' A„ ^ + A„ A„ ' + A„ ^ A n ' ^ dcj) n , 

where A^n denotes an amplitude with n external partons 
and I loops. d<f> n is the phase space measure for n par- 
tons. Taken separately, each of these contributions is di- 
vergent. Only the sum of all contributions is finite. To 




Fig. 2. Diagrams contributing to different primitive ampli- 
tudes. 

render the individual contributions finite, one adds and 
subtracts suitable pieces: 

{0) NNLO = 

J (<D n+2 da n % - O n +i o da [ ^l - O n o da^ 2) ) 

+ J (On+l da n % + On+l O dc^+l - O n O rf^ 1 ' 1 )) 

+ J (o n da n 2) + O rl o dan°' 2) + On o da n hl) ) . 

Here da n °+i is a subtraction term for single unresolved 
configurations of Born amplitudes. This term is already 
known from NLO calculations. The term dal°' 2 ' ) is a sub- 
traction term for double unresolved configurations. Fi- 
nally, dan is a subtraction term for single unresolved 
configurations involving one- loop amplitudes. 

To construct these terms the universal factorization 
properties of QCD amplitudes in unresolved limits are es- 
sential. QCD amplitudes factorize if they are decomposed 
into primitive amplitudes. Primitive amplitudes are de- 
fined by a fixed cyclic ordering of the QCD partons, a 
definite routing of the external fermion lines through the 
diagram and the particle content circulating in the loop. 
Fig. 2 shows three one-loop diagrams for e + e~ — > 3 jets 
contributing to different primitive amplitudes. One-loop 
amplitudes factorize in single unresolved limits as [24-30] 

^^Sing^-^x + SW 1 - 1 ).^. (1) 

Tree amplitudes factorize in the double unresolved limits 
as [31-37] 

4 0) = Sing(°> 2 > • 4°V 

To discuss the term dctn' 2 ^ let us consider as an example 
the Born leading-colour contributions to e + e~ — > qggq, 
which contribute to the NNLO corrections to e + e~ — > 
2 jets. The subtraction term has to match all double and 
single unresolved configurations. The double unresolved 
configurations are: 

- Two pairs of separately collinear particles, 

- Three particles collinear, 

- Two particles collinear and a third soft particle, 

- Two soft particles, 

- Coplanar degeneracy. 

The single unresolved configurations are: 

- Two collinear particles, 

- One soft particle. 
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Fig. 3. Splitting topology. 



It is convenient to construct dan^ as a sum over several 
pieces, 



topologies T 



Each piece is labelled by a splitting topology. An example 
is shown in fig. 3. The term (T) corresponding to 

the topology shown in fig. 3 approximates singularities in 
l/si2, l/(si2Si23) and part of the singularities in 1/s 2 



123- 



Care has to be taken to disentangle correctly overlapping 
singularities like l/(si2S23)- Details can be found in [14]. 



4 One-loop amplitudes with one unresolved 
parton 



Apart from dal?' 2 ^ also the term dan' 1 ', which approx- 
imates one-loop amplitudes with one unresolved parton, 
is needed at NNLO. If we recall the factorization formula 
(1), this requires as a new feature the approximation of 
the one-loop singular function Sing^ 1 ' 1 **. The correspond- 
ing subtraction term is proportional to the one- loop 1 — > 2 

splitting function V^' q) a _ fbc - An example is the leading- 
colour part for the splitting q — ► qg: 



(i,i) 



(1,0) q^qg,lc,corr 



iip(o.i) + s-i Cr 
q £ ' q->qg' TU £ w 

2 1 



ffl,corr(l/, Z) Vq-Xqg + h J*e [1 " - V )(l ~ z)] 

Sijk V 

This term depends on the correlations among the remain- 
ing hard partons. If only two hard partons are correlated, 
gi is given by 



gi,intr(y, z) 



1 



r(l + e)r{l-e) 



1 - z 



1 



-(1 - yfzf 2 F 1 (e, e, 1 + e; (1 - y)(l - z))} . 



For the integration of the subtraction terms over the un- 
resolved phase space all occuring integrals are reduced to 
standard integrals of the form 



l l 
dy y a (l - y) 1+c+d J dz z c (l - z) d [1 - z(l 
o o 
2 Fi (e,e; 1 + e; (1 - y)z) = 
r(i + a)r(l + d)r{2 + a + d + e)r(l + e) 
r(2 + a + d)r(e)r(e) 



y)Y 



E 



£ti + e)r{j + e)r(j + 1 + c) 
j=Z r ^ + ^ r ti + 1 + £ ^ r{ 3 + 3 + a + c + d + e)" 

The result is proportional to a hyper-geometric functions 
4^3 with unit argument and can be expanded into a Lau- 
rent series in e with the techniques of [38,39]. For the 
example discussed above one finds after integration [15]: 



V, 
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where L — ln(sjjfe//x ). The parameter p specifies the vari- 
ant of dimensional regularization: p = 1 in the conven- 
tional or 't Hooft-Veltman scheme and p = in a four- 
dimensional scheme. 



5 Outlook 

In this talk I reported on the subtraction method to can- 
cel infrared divergences at NNLO. The set-up involves two 
new types of subtraction terms, da$' 2 ^ and ddn' 1 ^ ■ The 
former approximates double unresolved configurations of 
tree amplitudes with n + 2 partons, whereas the latter ap- 
proximates one-loop amplitudes in single unresolved lim- 
its. Decomposing the QCD amplitudes into partial and 
primitive amplitudes, the appropriate subtraction terms 
have been constructed. Furthermore, the analytic integra- 
tion over the unresolved phase space has been performed 

for all terms contributing to dotn' 1 ^. Once the correspond- 
ing analytic integration has been done for da4°' 2 ^ the sub- 
traction method at NNLO is complete and can used for 
fully differential programs at NNLO. 
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